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ABSTRACT 

We present medium resolution (A/ AA « 2500) optical spectroscopy of SN 1999by in NGC 2841 
made around its light maximum. The depth ratio of the two Si II features at 5800 A and 6150 
A being i?(SiII)w 0.63 at maximum indicates that this SN belongs to the peculiar, sub-luminous 
SNe la. Radial velocities inferred from the minimum of the 6150 A trough reveal a steeper decline 
of the velocity curve than expected for "normal" SNe la, consistent with the behavior of published 
VRI light curves. A revised absolute magnitude of SN 1999by and distance to its host galaxy 
NGC 2841 is estimated based on the Multi-Color Light Curve Shape (MLCS) method, resulting 
in My (max)=— 18.06 ±0.1 mag and d — 17.1 ± 1.2 Mpc, respectively. An approximative linear 
dependence of the luminosity parameter A on i?(SiII) is presented. 



Subject headings: stars: supernovae: general — 
distances 

1. Introduction 

In the last decade it became evident that some 
of the observed Type la supernovae are intrinsi- 
cally sub-luminous compared to the majority of 
so-called "normal" SNe la. SN 1991bg (Fihppenko 
et al., 1992) is now considered as their prototype, 
and the sample of well-observed events includes 



^ Based on observations obtained at David Dunlap Ob- 
servatory, University of Toronto, Canada 

■^Visiting Astronomer, David Dunlap Observatory, Uni- 
versity of Toronto, Canada 



stars: supernovae: individual (SN 1999by) — stars: 

SN 1992K (Hamuy et al., 1994), SN 1997cn (Tu- 
ratto et al., 1998) and SN 1998de (Modjaz et al., 
2000). Beside being about 1 magnitude fainter 
than normal SNe la at maximum, sub-luminous 
SNe la also exhibit peculiar colors and spectra. 
For example, the intrinsic (unreddened) {B — T^)o 
color around maximum is 0.7 — 0.8 magnitude red- 
der than normal SNe la, reaching [B — F)o — 1.5 
mag at 10 days after maximum light. The spec- 
troscopic peculiarity is represented by i) the ap- 
pearance of strong Ti II band in the blue region 
and n) the increase of the ratio of the depths 
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SN 1999by (NGC 2841) 




Fig. 1. — Left: spectra of SN 1999by, shifted vertically for better visibility. Epochs are indicated on the 
right side of each spectrum. Right: CCD-image of SN 1999by (marked) taken with the Im RCC-telescope 
at Piszkesteto station of the Konkoly Observatory on May 21, 1999 (north is up and east is to the left). 



of Si II absorption troughs at 5800 A and 6150 
A (i?(SiII)) (e.g. Nugent et al., 1995). More- 
over, sub- luminous SNe la have steeply declin- 
ing light curves. Their decline rate parameter is 
Ami5{B) > 1.9 instead of 1.0 — 1.5 as for normal 
SNe la (e.g. Phillips et al., 1999). Because both 
i?,(SiII) and Ami5{B) correlate with the peak lu- 
minosity, the observational data indicate that sub- 
luminous SNe la have less energetic explosions and 
less ejected Ni masses than normal SNe la, al- 
though the correlation between the different quan- 
tities cannot be described with a single parameter 
(Hamuy et al., 1996, Hatano et al., 2000). 

In this paper we present medium resolution 
spectra of the sub-luminous la-type supernova 
SN 1999by occured in NGC 2841. Chronicles of 
the discovery and the earliest observations can be 
found in lAU Circulars No. 7156-7159. Gerardy 
& Fesen (1999) pointed out very soon that the 
spectrum of SN 1999by is of Type la, while Gar- 
navich et al. (1999) noted that SN 1999by belongs 
to the group of sub-luminous SNe la based on the 
i?(SiII) line depth ratio. Early-phase spectra of 
such SNe are available only for SN 1991bg and 
SN 1997cn. Thus, the spectroscopic observations 
made around maximum light may supply interest- 
ing details about peculiar SNe la. Very recently 



Howell et al. (2001) presented spcctropolarimot- 
ric evidence for the strong (« 20%) asphericity of 
the ejecta of SN 1999by, suggesting a connection 
between the observed asymmetry and the mecha- 
nism producing sub-luminous SNe la. 

Our new observations are described in Sect. 2, 
while their interpretation is given in Sect. 3. Sect. 4 
summarizes the results of this paper. 

2. Spectroscopic observations 

SN 1999by was observed with the 1.88 m tele- 
scope and Cassegrain spectrograph of the David 
Dunlap Observatory on three nights between May 
10 and 21, 1999. The 150A grating in 2nd order 
with an order separation filter in the stellar beam 
was applied, resulting in a reciprocal dispersion of 
1.22 A/pixel on the CCD detector. The recorded 
spectra cover about 1250 A centered on 6000 A. 

Reduction was performed by standard routines 
in IRAF^. The wavelength scale was calibrated 
with FeAr spectral lamp exposures. The back- 
ground light due to the host galaxy was approxi- 

^IRAF is distributed by NOAO which is operated by 
the Association of Universities for Research in Astronomy 
(AURA) Inc. under cooperative agreement with the Na- 
tional Science Foundation 
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mated with a quadratic polynomial perpendicular 
to the dispersion axis and subtracted from the SN 
aperture at each wavelength. 

Because of the rarity of photometric nights at 
the DDO site (in the proximity of Toronto), pre- 
cise flux calibration of the spectra was not pos- 
sible. However, an approximate correc;tion for 
the CCD spectral response has been determined 
by obtaining the spectra of three nearby field 
stars (#2, #3, #4 of Skiff, 1999; sec Fig.l of 
Toth & Szabo, 2000). Their {B - V) indices 
(0.59, 0.75, 0.60 mag, respectively) corresponds to 
Te// w 5500 - 6000 K according to the tables 
of Schmidt-Kaler (e.g. Carroll & Ostlie, 1996). 
The measured continua of these stars could be ap- 
proximated by a smooth polynomial. The real 
continuum fluxes (assuming log 5 = 4) were es- 
timated from model spectra of Buser fc Kurucz 
(1992). Then, the wavelength-dependent sensitiv- 
ity of the used CCD spectrograph was calculated 
by dividing the measured continua with the cho- 
sen model continuum flux values. The resulting 
function turned out to be smooth and slightly de- 
creasing to the red. We have looked for any sudden 
drop in the sensitivity function that might distort 
the depths of spectral lines and found nothing. 

The SN spectra were then corrected by divid- 
ing the measured spectra with the approximate 
sensitivity function described above. Because the 
used comparison stars were not spectroscopic stan- 
dards, the corrected "continuum" slopes of the SN 
spectra may still be slightly distorted. Neverthe- 
less, the SN spectra calibrated in such a way are 
acceptable for a comparative analysis described in 
Sect.3.1. 

The reduced spectra are collected and plotted 
in Fig.l. The epochs of the observations are in- 
dicated on the right side of each spectrum. The 
phases relative to B-maximum were determined by 
adopting the result of Bonanos, et al. (1999) that 
the maximum in B occured on JD 2451308.5 (May 
9/10, 1999), which was also confirmed by photom- 
etry of Toth & Szabo (2000) (note that the JD of 
maximum light was misprinted in both papers). 

Although the spectral coverage is rather small, 
the appearance of the strong Si II line around 6150 
A unambigously identifies a Type la SN. The other 
strong feature, the absorption trough at 5800 A is 
also attributed to Si II. The relative strength of 
this line to that of the 6150 A fine i?(SiII) (Nu- 



gent et al., 1995) suggests that SN 1999by is a 
sub- luminous event, as mentioned in Sect.l. This 
is illustrated in Fig. 2, where the maximum-light 
spectrum of SN 1999by is plotted together with a 
pre-maximum spectrum of SN 1998aq (Vinko et 
al., 1999) obtained with the same instrument and 
setup (unfortunately, no spectrum of SN 1998aq 
closer to its maximum light was available for us). 
The spectra in Fig. 2 were normalized to the local 
"continuum" for illustrative purpose. It is visible 
that the 5800 A line is much weaker in SN 1998aq 
than in SN 1999by. This line gets stronger in pre- 
maximum spectra of SNe la when approaching the 
maximum, however, in "normal" SNe the increase 
of i?(SiII) is only 0.04-0.06 (Riess et al., 1998b). 
Thus, the 5800 A feature was not much stronger 
in the spectrum of SN 1998aq at maximum. 

3. Data analysis and interpretation 
3.1. Line depths and radial velocities 

We have calculated the value of i?(SiII) of 
SN 1999by following the definition of Nugent et 
al. (1995). i?(SiIl) = 0.63 was obtained for the 
May 11 spectrum, while R = 0.66 was determined 
for the May 14 one, in accord with the result of 
Garnavich et al. (1999) (see also Hatano et al., 
2000) . This is more than a factor of 2 larger than 
its usual value oi R = 0.25 ±0.1 for normal SNe 
la (Riess et al, 1998b, Hatano et al., 2000). 

Very recently Howell et al. (2001) published 
three spectra of SN 1999by made just before 
our observations. From their Fig.l we estimated 
i?(SiII) = 0.68, 0.70, 0.66 at epochs r = -2, -1, 0, 
respectively. These agree very well (within a few 
hundredths) of the results derived from our spec- 
tra, suggesting that our observations probably do 



Table 1: Radial velocities of SN 1999by from the 

minimum of the Si II line. The velocity of NGC 
2841 {vgx = 638 kms~^) has been added to get 
rest-frame radial velocities of the SN. 



JD 


T 


Vr 


(-2400000) 


(day) 


(kms"""^) 


51309.6 


+1 


10174 


51312.5 


+4 


9560 


51319.6 


+11 


8946 
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Fig. 2. — Comparison of spectra of SN 1999by 
(bottom) and SN 1998aq (top) around maximum 
light. The latter spectrum was shifted by 0.5. The 
relative strength of the 5800 A Si II trough in SN 
1999by is apparent. 
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Fig. 3. — Radial velocity curve of SN 1999by com- 
puted from the minimum of the Si II 6150 A line 
(boxes), compared with those of SN 1998aq (dia- 
monds) and SN 1994D (pluses). See text for ref- 
erences. 



not contain wavelength-dependent systematic er- 
rors distorting the line-depth ratio. 

i?(SiII) correlates with the luminosity param- 
eter A introduced by Riess et al. (1996) in the 
Multi-Color Light Curve Shape (MLCS) method. 
A calibration of the linear dependence of i?(SiII) 
on A was given by Riess et al. (1998b), re- 
sulting in A = 2.33i? - 0.55(±0.22) for -|-1 day 
phase, while for +4 days the slope and zero point 
changes to 2.06 and —0.46, respectively. Substi- 
tuting the i?(SiII) values given above one can get 
A = 0.90 ± 0.23 for the luminosity parameter of 
SN 1999by (the given error takes into account the 
quoted uncertainty of the i? - A calibration). For 
"normal" SNe la A » 0, therefore, this high value 
suggests that SN 1999by is one of the most sub- 
luminous SNe la (Bonanos et al., 1999). Note that 
this spectroscopic A is not suitable for describing 
the multi-color light curve of SN 1999by, which 
requires A > 1.0. The light curve decline parame- 
ter. Amis (5) = 1-87, is also indicative of the low 
luminosity (Bonanos et al., 1999). 

Radial velocities were inferred from the absorp- 
tion minima of the 6150 A Si II trough (fitted by 
a low-order polynomial around minimum). Table 
1 lists the resulted radial velocities, corrected for 
the host galaxy center recession velocity Vgx = 638 
kms~^ (as given in the NASA/IPAC Extragalac- 
tic Database). As the anonymous referee pointed 
out, SN 1999by was closer to the edge of NGC 2841 
and could have significant orbital velocity in ad- 
dition to the center-of-galaxy radial velocity. The 
rotation velocity of NGC 2841 at distances greater 
than 10 kpc from its center is « 300 kms^^ (Sofue 
et al., 1999). Thus, the uncertainty of the radial 
velocities is at least ± 300 kms^^ due to the possi- 
ble orbital velocity as well as the large width of the 
Si II line. Nevertheless, this uncertainty is only a 
few percent of the Si II expansion velocity. 

Fig. 3 compares the radial velocities of SN 1999by 

with those of SN 1998aq (Vinko et al.,1999) and 
SN 1994D (Patat et al.,1996). The decline of the 
radial velocities of SN 1999by is much stronger 
than in the other two cases. SN 1991bg showed 
similar behavior (Filippenko et al.,1992), while the 
other sub-luminous SNe la were not well covered 
spectroscopically around maximum. The stronger 
velocity decline is also in good agreement with the 
recent result of Hatano et al. (2000) that sub- 
luminous events have lower velocities at 10 days 
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after maximum than do normal SNc la. This may 
support the hypothesis that sub-luminous SNe la 
arise from sub-Chandrasekhar mass progenitors, 
thus have less energetic explosion than normal 
Type la events (see Modjaz et al., 2000 for recent 
review of explosion models). 

3.2. Absolute magnitude and distance 




290 300 310 320 330 340 350 360 370 
JD -2451000 

Fig. 4.— MLCS-fit to published BVRI data of 
SN 1999by. 



Application of the Multi-Color Light Curve 
Shape (MLCS) method for SN IQQQby has been 
published very recently by Toth & Szabo (2000). 
They used only the V light curve for template 
fitting and {B — V) data for the reddening de- 
termination, and revealed A > 1.6 mag. This 
may indicate that SN 1999by is extremely sub- 
luminous. For example, SN 1991bg, the prototype, 
had A = 1.44±0.1 (Riess, Press & Kirshner, 1996) 
which is known as one of the least luminous SNe 
la to date. 

We re-analysed the BVRI photometric dataset 
published by Toth & Szabo (2000) supplemented 
by the data of Hanzl (1999) using the original 
MLCS template curves of Riess, Press & Kirsh- 
ner (1996). Because the photometric data are in- 



homogencous, systematic errors may distort the 
parameters derived from them. Thus, the results 
below can be considered only as crude estimates. 

Although Riess et al. (1998a) developed a sec- 
ond version of MLCS (MLCS-2, based on an ex- 
tended list of calibrating SNc), the updated tem- 
plate vectors have not been published yet. How- 
ever, recently Vinko et al. (2001) demonstrated 
that for nearby SNe, even with moderate redden- 
ing, the MLCS-1 and MLCS-2 methods produce 
consistent results, if the maximum magnitude of 
the fiducial (A — 0) light curve is set properly. 
We adopted M{?"^ = -19.46 mag for the fiducial 
V light curve (Riess, personal communication), 
while the original MLCS-1 method was based on 
]\/[rno.x _ _X9_36 mag. The difference has been 
taken into account by adding -1-0.1 mag to the dis- 
tance modulus computed with MLCS-1. 

For SN 1999by, E{B -V) = 0.01 was adopted, 

consistent with the all-sky reddening map of 
Schlegel et al. (1998) (Toth & Szabo derived 
E{B-V) = 0.05±0.02 from their SN photometry). 
The template vectors were fitted simultaneously 
to all BVRI data, i.e. the residual of all data were 
combined in a single x^-function. The fitting was 
restricted for r > —2 days, because the available 
pre-maximum data of SN 1999by could not be ad- 
equately fitted together with the post-maximum 
data (see also Toth & Szabo, 2000). Also, the min- 
imum photometric error was increased to ±0.15 
mag taking into account some possible system- 
atic uncertainties in the used dataset (due to e.g. 
errors in the standard transformation). 

The simultaneous fitting of the template curves 
resulted in A = 1.4±0.1 and /io = 31.16±0.15 mag 
for the luminosity parameter and the reddening- 
free distance modulus, respectively. The given un- 
certainty of the distance modulus is only the rms 
error of the fitting, disregarding the possible sys- 
tematic errors in the absolute magnitudes of the 
template light curves. The quality of the fit can 
be judged from Fig. 4. It is seen that the fitting is 
not perfect, however, we believe that this is close 
the optimal result that can be achieved from the 
present datasets. The minimum of the reduced 
chi-squared that could be reached was ~ 
Still, the solution may be systematically incorrect, 
because for such under- luminous SN the template 
vectors in MLCS-1 are based only on the light 
curves of SN 1991bg, and the light variation of 
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R(Sill) 



Fig. 5. — Dependence of the light curve parameter 

A on i?(SiII) using data of well-obscrvcd SNc la 
(see Table 2 for references). The dotted line rep- 
resents the relation given by Riess et al. (1998b), 
while the solid line shows our new calibration (see 
text). 

SN 1999by may differ substantially from that. The 
large deviation of the pro-maximum points may be 
due to such inconsistence. Also, photometric er- 
rors higher than we adopted (±0.15 mag) cannot 
be ruled out. 

The photometric A is significantly higher than 
the value found from spectroscopy (A = 0.90, 
Sect. 3). It was not possible to reach an adequate 
fit with A = 0.9 given by spectroscopy. Thus, the 
current spectroscopic determination of A does not 
give a suitable result for SN 1999by, which is prob- 
ably true for other extreme sub-luminous SNe la 
as well. An improved i?(SiII) — A calibration is 
presented in the next section. 

Although sub-luminous SNe la are usually elim- 
inated from the sample of "good" distance indi- 
cators, the MLCS method gives information on 
their distances as well. The new distance mod- 
ulus of SN 1999by given above corresponds to 
17.1 ±1.2 Mpc geometric distance, which is sig- 
nificantly higher than d = 10 — 12 Mpc given by 
Toth & Szabo (2000) and TuUy (1988). Note that 



the template vectors of the MLCS-1 method were 
calibrated from a training set of SNe adopting 
Tully-Fisher- and surface-brightness-fluctuation 
distances, but the final absolute magnitudes were 
shifted to match the Cepheid-based distance scale 
(Riess, Press & Kirshner, 1996). Thus, the higher 
distance given above simply reflects the differ- 
ence between various distance scales, and may 
still contain a systematic error. Both Cepheid- 
and Tully-Fisher distances undergo signiflcant re- 
visions nowadays (e.g. Gibson et al., 2000, Sakai 
et al., 2000), so this question can be answered only 
when the problems of the different distance scales 
are solved. 

3.3. Comparison with other SNe la 

The absolute magnitudes of SNe la correlate 
with i?(Sill) as was presented very clearly by Nu- 
gent et al. (1995). Although Hatano et al. (2000) 
demonstrated that this correlation is not a one- 
parameter sequence, the dependence of M™ ''^ on 
i?(SiII) may still be a somewhat useful relation 
for predicting SNc luminosities from spectra ob- 
served around maximum (at least as a first ap- 
proximation). Now the correlation can be slightly 
improved using the new data of sub-luminous SNe 
la. 

Table 2: ii(Sill) and A parameters of SNe la at 
maximum. 



SN 


i?,(SiII) 


A 




Ref. 


1994ae 


0.12 


-0.13 


0.05 


a,b 


1991T 


0.14 


-0.56 


0.31 


a,b 


1981B 


0.16 


0.00 


0.32 


a,b 


1990N 


0.16 


-0.28 


0.31 


a,b 


1998bu 


0.23 


0.02 


0.05 


d,e 


1989B 


0.29 


0.04 


0.31 


a,b 


1992A 


0.38 


0.44 


0.13 


a,b 


1986G 


0.53 


0.41 


0.31 


a,b 


1991bg 


0.62 


1.44 


0.10 


a,b 


1997cn 


0.63 


1.50 


0.05 


c 


1999by 


0.63 


1.40 


0.10 


PP 



References: a) Nugent et al. (1995); b) Riess, 
Press & Kirshner (1996); c) Turatto et al. (1998); 
d) Jha et al. (2000); e) Hatano et al. (2000); pp) 
present paper 

We have collected values of the light curve 
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parameter A that were derived via the MLCS 
method and the values of i?(SiII), both at V^— maxi- 
mum. These are Hsted in Table 2. Because A 
measures the deviation of the V light curve of a 
particular SN from the fiducial SN la light curve 
at maximum, it is also directly proportional to the 
maximum absolute magnitude M™ However, 
at present, A seems to be a more useful parameter 
than M™"^, because it is determined directly from 
the shape of the (multi-color) light curves, there- 
fore it is not affected by the adopted zero-point of 
the distance scale. 

The result is plotted in Fig. 5, where the contin- 
uous line represents the approximate relation 

A = 3.5 i?(SiII) - 0.8 (1) 
±0.44 ±0.17 

while the dashed line shows the relation given by 
Riess et al. (1998b) (see Sect. 3.1). It is seen that 

the Riess et al. calibration gives systematically 
lower A for the most sub-luminous SNe la. The 
new relation fits these SNe better, thus it may 
be used to predict A from the measurement of 
i?(SiII) around maximum. Such prediction may 
be useful, for example, in the "snapshot" distance 
estimate method (Riess et al., 1998b) when only a 
few nights' data are available for a particular SN. 

From Fig. 5, it may also be interesting that the 
most sub-luminous SN 1991bg, SN 1997cn and 
SN 1999by represent a small group with very con- 
sistent line-depth ratio and A (or absolute mag- 
nitude). It may mean that sub-luminous SNe la 
might be better distance indicators than previ- 
ously thought, if their sub-luminous status is prop- 
erly identified from spectroscopy and their differ- 
ence from normal SNe la is taken into account. 
It is known (e.g. Hamuy et al., 1994) that such 
sub-luminous SNe la may be more frequent events 
than their discovery rate suggests, thus, it is pos- 
sible that many more such objects will be detected 
in the extensive supernova-search programs. How- 
ever, the present sample is obviously too small to 
draw a certain conclusion. 

4. Summary 

The results of this paper can be summarized as 
follows: 

Three spectra of the sub-luminous Type la 
SN 1999by around maximum were obtained at 



DDO, centered on the Si II trough at 6150 A. The 
i?(SiII) line depth ratio (Nugent et al., 1995) is 
determined to be 0.63 at maximum, indicative of 
an intrinsically sub-luminous SNe la. 

Radial velocities of the 6150 A Si II line were 
determined from measuring the Doppler-shift of 
the line core. It is shown that SN 1999by exhibited 
stronger decline of its radial velocity curve than 
other, spectroscopically normal SNe la. This may 
support the low-energy, sub-Chandrasekhar mass 
explosion proposed for such SNe. 

Absolute magnitude and distance of SN 1999by 
have been derived via the MLCS method. The as- 
sumption that the reddening is E{B — V) = 0.01 
led to M{?''^=-18.06±0.1 mag and d = 17.1 ± 1.2 
Mpc. The maximum absolute magnitude as well 
as the line depth ratio i?(SiII) of SN 1999by are in 
very good agreement with those of other extreme 
sub-luminous SNe la, SN 1991bg and SN 1997cn. 
An approximative linear relation has been deter- 
mined to predict the light curve parameter A of 
SNe la over a wide range of observed i?(SiII). 
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